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Websites offering Free Online Courses: 


Learn. Think. Do. 

Advance Your Education with Free College Courses Online: www.udacity.com 

Free Online Courses: www.coursera.org 

Explore Free Courses from edX Universities: www.edX.org 

Feam almost anything for free: www.khanacademy.org 

Soziale Plattform fur Interaktive Online-Kurse zur Informationstechnologie (deutsch) 
www.openhpi.de 

Download die Vortrage der Lehrertage an der Universitat Wien fur viele Jahre: 
http://www.oemg.ac.at/DK/Didaktikhefte/index.html 


The DERIVE-NEWS archive: 

https://www.jiscmail.ac.uk/cgi-bin/webadmin? AO=derive-news 


Information from Djordje Kadijevich, Serbia: 

Dear Colleagues, 

Hope you are fine. 

Just to tell you about my new book published by Routledge: 

Improving Computer Science Education, suitable for courses 
on didactics of informatics (didactics of computer science). 

More information about this book can, for example, be found at the Publisher web site at 

www . routledge . com/books/ details/ 97 804 15 64537 9/ 

For on-line orders at this site, you can enter code ERJ60 to receive a 20% discount. 
Could you please circulate this information among interested colleagues and students. 
Kind regards, 

Djordje 
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Liebe DUG-Mitglieder, 

DNL#89 hat einige Zeit zur Fertigstel- 
lung benotigt. Er war sehr arbeits- und 
zeitaufwandig, da die meisten Beitrage 
erst ins Englische ubersetzt werden 
mussten. Herzlichen Dank an Herrn Erik 
van Lantschoot, der nicht nur den schonen 
Beitrag zu einem historischen Bauwerk 
seines Geburtsorts sondern auch die 
Ubersetzung beisteuerte. Fur die anderen 
Ubersetzungen bin ich ganz alleine ver- 
antwortlich. 

Ich musste einige - sehr interessante - 
Beitrage zum User Forum fur den nachs- 
ten DNL aufheben. Alle Anfragen haben 
ausfuhrliche Beantwortungen nach sich 
gezogen. Dies wiirde den Rahmen des DNL 
bei weitem sprengen. 

So hat mich die Frage von Herrn Korner 
nach dem „Wert" der Polynomdivision sehr 
lange und ausfuhrlich beschaftigt. Hof- 
fentlich konnen Sie Teile meiner Antwort 
auch gebrauchen. Wenn Sie noch weitere 
Kommentare beitragen wollen, sind Sie 
herzlich dazu eingeladen. 

Kiirzlich kam die Meldung, dass die 
DERIVE-NEWS list eingestellt wird, wenn 
sich kein neuer Administrator melden 
wiirde. Gliicklicherweise hat sich Lars 
Erup (Kanada) dazu bereit erklart. Herzli- 
chen Dank dafiir und viel Erfolg! Auf Sei- 
te 47 finden Sie seine Nachricht und die 
Internetverbindung zu dieser News 
Group. Bitte nut zen Sie diese weitere 
Kommunikationsmdglichkeit. 

Am Ende dieses DNL finden Sie den Fol- 
der zur TIME 2014, die im schonen 
Krems, Niederosterreich, stattfinden 
wird. 

Ich verbleibe mit den besten GriiBen bis 
zum ntichsten Mai 

Ihr 


Dear DUG Members, 

DNL#89 needed some time to be fin- 
ished. Most of the contributions are 
translations from German. Many thanks to 
Erik van Lantschoot. He not only contrib- 
uted his fine article describing a detail of 
a historic landmark of his birthplace but 
also the translation. I am full responsible 
for all other translations. 

I must leave some - very interesting - 
contributions for our User Forum for the 
next newsletter. They all require ex- 
tended answers which would be far be- 
yond the space of this DNL. 

Take Klaus Korner's question for the 
"value" of polynomial long division as an 
example. I hope that you can use at least 
parts of my answer. Additional comments 
would be very welcome. 

It was just recently when I received a 
mail that the DERIVE News list would be 
discontinued after many years of its exis- 
tence if nobody from within the news 
group would volunteer as list administra- 
tor. Fortunately Lars Erup (Canada) took 
over the job. Many thanks for this and 
much success in the future! You can find 
his message on page 47 together with the 
internet connection to the DERIVE News 
Group. Please use this additional possibil- 
ity for communication within the DERIVE- 
CAS community. 

At the end of this DNL you can find the 
information folder for TIME 2014 which 
will be held in July 2014 in lovely Krems, 
Lower Austria. 

I remain with my best regards 

until DNL#90 

Yours 



Download all ZWL-Derive- and Tl-flles from 

http : / / www . austromath . at/ dug/ 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 77 User Group. 
It is published at least four times a year 
with a contents of 40 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE, 
77-CAS and other CAS as well to create a 
group to discuss the possibilities of new 
methodical and didactical manners in 
teaching mathematics. 


Editor: Mag. Josef Bdhm 
D'Lust 1, A-3042 Wurmla, Austria 
Phone: ++43-(0)660 3 1 36 365 

e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 

Next issue: June 2013 


Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 
Wonderful World of Pedal Curves, J. Bohm, AUT 
Tools for 3D-Problems, P. Luke-Rosendahl, GER 
Hill-Encryption, J. Bohm, AUT 

Simulating a Graphing Calculator in DERIVE, J. Bohm, AUT 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 

Graphics World, Currency Change, P. Charland, CAN 

Cubics, Quartics - Interesting features, T. Koller & J. Bohm, AUT 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ / Pierre Charland 's Graphics Gallery 

BooleanPlots.mth, P. Schofield, UK 

Old traditional examples for a CAS - what’s new? J. Bohm, AUT 

Truth Tables on the Tl, M. R. Phillips, USA 

Where oh Where is It? (GPS with CAS), C. & P. Leinbach, USA 

Embroidery Patterns, H. Ludwig, GER 

Mandelbrot and Newton with DERIVE, Roman Hasek, CZK 

Some Projects with Students, R. Schroder, GER 

Dirac Algebra, Clifford Algebra, D. R. Lunsford, USA 

Laplace Transforms, ODEs and CAS, G. Picard & Ch. Trottier, CAN 

A New Approach to Taylor Series, D. Oertel, GER 

Cesar Multiplication, G. Schodl, AUT 

Henon & Co; Find your very own Strange Attractor, J. Bohm, AUT 
Rational Hooks, J. Lechner, AUT 

Simulation of Dynamic Systems with various Tools, J. Bohm, AUT 
An APL-like SHAPE function in DERIVE 6, D. R. Lunsford, USA 
Recursive Series of Numbers, Polynomials and Functions, D. Halprin, AUS 
Inspirations from other sources, J. Bohm, AUT 
Teaching Calculus, G. Savard, CAN 
and others 


Impressum: 
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Richtung: Fachzeitschrift 
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A Mathematical Model for Snail Shells (4) 

Piotr Trebisz, Germany 


Der Ansatz, wonach man das Selbstdurchdringungsproblem durch Verwendung krummli- 
niger Koordi natensysteme analytisch losbar macht, kann naturlich auch auf koni- 
sche Schneckenhauser ubertragen werden. Konische Schneckenhauser dieser Art nen- 
ne ich „Trebi sz-Spi ral en - Typ 4". 

Es handelt sich dabei aber um mehr als nur um einen mathemati schen Trick oder 
reine Spielerei! Fakt ist dass die sich so ergebenden Schneckenhauser insgesamt 
schoner und eleganter sind und - mehr noch als ihre Gegenstucke vom Typ 2 - aus 
einem Guss erscheinen, weil bei ihnen tatsachlich samtliche geometri schen Eigen- 
schaften durch die logarithmische Spirale bestimmt werden. 

Im final en 4. Teil meiner Bei tragsrei he werde ich zeigen wie das korrespondi e- 
rende krummlinige Koordi natensystem konischer Schneckenhauser aus der Spiralen- 
Glei chung hergeleitet wird. Die Vorgehenswei se ist dabei analog zu Trebi sz- 
Spi ral en von Typ 3 . 

Ausgangspunkt unserer Uberlegungen ist die Gleichung einer konischen logarith- 
mi schen Spirale mit den 3 Parametern b (Basis), t (Lange) und m (Proportiona- 
litatsfaktor) und ihre Tangente. 

Gesucht ist nun eine Anti -Spirale mit der Basis a, die die normal e Spirale in 
den Schni ttpunkten immer senkrecht schneidet. Ihr Proportionalitatsfaktor m ist 
identisch, so ist gewahrlei stet dass sie sich auf derselben Kegeloberflache win- 
det. (Siehe auch meinen Kommentar auf Seite 11, Josef) 


Um Platz zu sparen sind Piotrs Kommentare nur in (meiner) englischen Fassung, Josef. 


The concept solving the interpenetrating problem by using curvilinear coordinate systems can be 
transferred to conical snail shells, of course. I will call these snail shells “Type 4 Trebisz-Spirals ” . 

This is more than a mathematical trick or pure playing around. Fact is that the resulting snail shells are 
much more beautiful and elegant. They appear all of a piece - even more than the respective 
type 2 models do because all their geometric attributes are determined by the logarithmic spiral. 

In this 4 th - and final - contribution I will show how the corresponding curvilinear coordinate system 
of conical snail shells can be derived from the equation of the spiral. The procedure is the same as 
demonstrated for the Type 3 Trebisz Spirals. 

Important settings for the state variables are: 

[Precision := Exact, PrecisionDigits := 10] 

[Notation := Rational, NotationDigits := 10] 

[Branch := Principal, Exponential := Auto, Logarithm := Auto] 

[Trigonometry := Auto, Trigpower := Auto, Angle := Radian] 

[InputMode := Word, CaseMode := Sensitive] 

[b :c Real (1, oo) , m :e Real (0, oo) , t :e Real [0, oo)] 

[0 :e Real, :e Real [-7T, +7r], a :e Real (0, 1), x :e Real [0, »)] 

[Ri :e Real (0, oo) , R 2 :e Real (0, oo) , k :e Integer] 
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Starting point for our considerations is the equation of a conical logarithmic spiral defined by three 
parameters: b (base), t (length), and m (proportionality factor). An anti-spiral (base a) is needed which 
intersects the given spiral under 90°. Its prop factor is the same ( m ) which makes sure that both spirals 
are winding on the same conic surface. They intersect for t = 1 . 

SPIRALED, m, t) := [t, m. t- C0S(2 .jr.L0G(t, b)), m.t.SIN(2.n-.L0G(t, b))] 

d 

SPIRAL E_TANGENTE(b , m, t) := — SPIRALED, m, t) 

dt 

ANTI_S P I R AL E ( a , m, t) := [t, m- 1 ■ C0S(2 -ir- L0G(t , a)), m-t.SIN(2.ir.L0G(t, a))] 

d 

ANTI_S P I R AL E_TAN GE NTE ( a , m, t) := — ANTI_S P I R AL E ( a , m, t) 

dt 


SPIRALED, m, 1) 


[1. m, 0] 

. ANTI_SPIRALE(a, m, 1) . 


i 

i 1 

3 

0 

1 1 


Base a of the anti-spiral must be such chosen that both spirals intersect under 90°. This condition is 
given in the next equation which is then solved for a {a = BASIS (b , m) ). 

SPIRAL E_TANGENTE(b, m, 1) . ANTI_SPIRALE_TANGENTE(a , m, 1) = 0 

BASIS (b , m) := (SOLUTIONS(SPIRALE_TANGENTE(b , m, 1) . 


ANTI_SPIRALE_TANGENTE(a , m, 1) = 0, a, Real)) 

1 


2 2 2 

- 4-n -m / ( (m + l)*LN(b)) 

BASIS (b, m) = e 

We still miss the local y-axis. It must be perpendicular to the tangent of the anti spiral and in the inter- 
section points of the spirals perpendicular to the tangent of the given spiral, too. For the type 3 shells is 
this the normal of the plane which contains the planar spiral. Here the normal of the cone surface turns 
out to be an appropriate local r-axis. This normal line is constructed by choosing to every point t on 
the anti-spiral a corresponding point x on the x-axis such that the line connecting these two points and 
the tangent of the anti-spiral are perpendicular. The length of the connecting segment depends linearly 
on t. 

LINIE(a , m, t, x) := ANTI_SPIRALE(a, m, t) - [x, 0, 0] 

PUNKT(a , in, t) := (SOLUTIONS(LINIE(a , m, t, x) . ANTI_SPIRALE_TANGENTE(a , 


m, t) = 0, x, Real)) 

1 

ANTI_SPIRALE_ABSTAND(a , m, t) := LINIE(a , m, t, PUNKT(a, m, t)) 

2 (*) 

| ANTI_SPIRALE_ABSTAND(a , m, t)| = m-t-YOn +1) 
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We must get rid of this dependency if we will use this line as y-axis in our curvilinear system of coor- 
dinates. Now we have got much closer to our goal because we know both axes of the coordinate sys- 
tem. 


ANTI_SPIRALE_ABSTAND(a , m, t) 

ANTI_SPIRALE_NORMALE(a , m P t) := 

t 


X_ACHSE(b , in, t) := ANTI_SPIRALE( BASIS (b p m) , m p t) 

Y_ACHSE(h , m, t) := ANTI_S PI RALE_NORMALE( BASIS (b P m) , m p t) 

In these curvilinear coordinates we can define the surface ellipse of the shell with radius R \ in 
x-direction and a radius R 2 in y-direction. 

ELLIPSECb, m , R lp R 2 P <f) := X_ACHSE(b, m, 1 + Ri-COSOjO) + Y_ACHSE(b , m, 1 
+ Ri*COS(<j0) * R 2 *SIN(t))) 

The centre of the ellipse is at t = 1 . How to choose both axes? Ri should be chosen that the surface of 
the solid does not interpenetrate, and R 2 in such a way that the ellipse osculates the global x-axis 
(without intersecting it). 

Let’s concentrate first on R\ because it is calculated like the radius of a type 3 shell. The intersection 
points of both spirals are needed. They both are parameterized wrt to the length and they both have 
identical x-components. So we need all values for t which give identical y- and z-components. The 
solutions should be known from the type 3 shells. 

SOLUTIONS((SPIRALE(b, m, t)) = (X_ACHSE(b, m, t)) a SOLVE((SPIRALE(b , m, t)) 

2 2 3 

= (X_ACHSE(b, m, t)) , t, Real), t, Real) 

3 

222 2 22 222 2 
- 4«7T *m /((m + l)*LN(b) + 4*7T «m ) 4*7T «m /((m + l)«LN(b) + 

I L 0 , 1, b , b 

2 2 ‘ 

4«7T *m ) 

DERIVE returns only 4 special solutions but there are infinite many of them. By inspection it is possi- 
ble to give a general formula, k is an arbitrary integer. 

The difference of two intersections gives the “winding number”, i.e. the number of rotations after 
which the surface is osculating itself. 

2 2 2 2 2 2 
4*7t «m 'k/CCm + l)-LN(b) + 4«7 t «m ) 

SCHNITT(b , m, k) := b 

WINDUNG(b , m) :z LOG(SCHNITT(b , m, k + 1), b) - LOG(SCHNITT(b , m, k), b) 

2 2 
4*7t *m 

WINDUNG(b , m) := 

2 2 2 2 
(m + l)«LN(b) + 4-7T -m 


Now it is possible to calculate radius R\ in direction of the x-axis. 
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RADIUS JC(b, m) :z (SOLUTIONS(SCHNITT(b , til, k).(l + R-l) = 5CHNlTT(b , m, 


k + 1).(1 - Rx), Ri f Real)) 

1 

2 2 2 2 2 2 

4-tt -m /((in + l)-LN(b) + 4 -tt ) 

b - 1 

RADIUS^KCb, m) := * 

2 2 2 2 2 2 

4-7T -m /((m + l)-LN(b) + 4-/r -m ) 

b + 1 


It is not too difficult to find radius R 2 (in direction of they-axis). As the ellipse shall only osculate the 
global x-axis we need to know the distance of any arbitrary point on the anti-spiral to this axis along 
the local y-axis. The attentive reader will remember equation (*) from above (page 4). This distance 
has already been calculated. It scaled according to the scale of the local y-axis and we change this by 
the respective division: 

2 

| ANTI_SPIRALE_ABSTAND(a , m, t)| z + 1) 

| ANTI_SPIRALE^ABSTAND(a , m, t) | 

= t 

| ANTI_SPIRALE_NORMALE(a, m, t) | 

This means that the distance of any point of the anti-spiral from the global x-axis is exact t. Consider- 
ing this fact we can set up an equation in order to determine R 2 . This works as follows: 


1. We use a 2-dimensional system of coordinates which is a representative for the local curvilin- 
ear system with its x- andy-axis. As a consequence a diagonal withy = t corresponds with the global 
x-axis which shall be tangent to the ellipse. 


2. Analogous to the type 3 shell a circle with radius R\ is placed for t= 1. Is equation is: 

2 2 

y = ^(Ri - (t - i) ) 

But this is not the requested ellipse because it should have radius R 2 in y-direction. How to achieve 
this? This is not difficult at all: just scale the circle: 

2 2 R 2 

y = i/(Ri - (t - i) ) 

Ri 


This equation enables finding the intersection points of the ellipse and the diagonal y = t. 

r 2 2 

y(RADIUS_X(b, m) - (t - 1) ) 


RADIU5_Y_TEMP(b , m, R 2 ) := SOLUTIONS 


RADIUS_X(b , m) 


■R 2 = t, t, 




Real 


Then we choose R 2 in such a way that both intersection points coincide because then - and only then 
the global ellipse is osculating the x-axis (without intersecting it). 






P 8 


Piotr Trebisz: A Mathematical Model for Snail Shells 


( 4 ) 


D-N-L#89 


which can be further simplified to: 

r 2 2 2 2 2 2 ^ 
4*71 «m /((m + l)«LN(b) + 4«7i -m ) 


0 = 


^b 


- 1 


2 2 
VCCm + l)*LN(b) + 


2 2 2 2 2 2 
2 2 2 2*71 «m /((m + l)*LN(b) + 4«7i *m ) 

4«7i «m ) - 4«7T«m -b 


Unfortunately this equation cannot be solved for m. We have no other choice: we must apply again the 
well-tried Newton-Raphson-algorithm for determining m numerically. 

Function/(m,6) is the equation to be solved. g(m,b) is its first derivative, and h(b,m,n) is calculating 
the next iteration of the algorithm accurate to n decimal places. ( * } 


f(b, 


m) := 


r 2 2 2 2 
4‘7T *m /((m + 1) -LN(b) 


2 2 

+ 4«7T ‘fn ) 



2 

•7((m 


2 

+ l)-LN(b) 


2 2 

+ 4-7T »m ) - 


2 2 2 2 2 2 
2 2 *7T *m /((m + l)*LN(b) + 4-7T *m ) 

4*7T-m -b 


d 

g(b, m) := — f(b, m) 
dm 

h(b, m, n) := 

Prog 

6m :€ Real (0, oo) 

WRITE(APPEND("m := ", APPROX(STRING(m) , n))) 

5m := m - APPR0X(f(b, m) , 2*n)/APPR0X(g(b, m) , 2-n) 
RETURN APPROXCSm, n) 


For finding approximations of the m-zeros of/(/?,m) we need a suitable initial value. We will investi- 
gate the run of the function. 



We see that the function starts increasing and then after passing its maximum point turns to minus 
infinity. 


See our dialogue at the end of Piotr’s contribution. Josef 
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Can we give a general statement about its slope behavior which does not depend on its base bl 
Yes, we can! First of all ,J[b,m= 0) = 0, and the limit —> oo) = -oo. 

[ f (b, 0), lim f(b, m)l = [0, -~] 
m^oo J 


This information is useful for finding a suitable initial value mo for the approximation algorithm. 

It is easy! We start with m 0 = 1 and then increase m 0 until /(/j,» 7 0 ) becomes negative. Then we go 
backwards to the zero applying the Newton-Raphson-algorithm. An implementation of this algorithm 
could look as follows: 

M(b, n) := 

Prog 

m 0 :e Real (0, oo) 
fo :e Real 

V :e Vector 
m 0 := 1 
Loop 

f 0 := APPROX(f(b, m 0 ), 2-n) 

If f 0 > 0 
m 0 := m 0 + 1 
If f 0 < 0 
exit 

RETURN m 0 

V := ITERATES(h(b, m, n) , m, m 0 ) 

RETURN ViDIM(V) 


Now we are really done and our award is the presentation of wonderful conical snail shells defined by 
three parameterizations. 

SchneckenHausX_U(b, n, 0, tp) := SchneckenTorusX_U(b, M(b, n), 0, tp) 

SchneckenHaus>L_W(b, n, 0, tp) := SchneckenTorusX_W(b, M(b, n), 0, tp) 

SchneckenHausX_L(b, n, t, tp) := SchneckenTorusX_L(b, M(b, n) , t, tp) 

(U for “rotation number” (= Umdrehungszahl), W for “winding number” (= Windungszahl) 
and L for “Length” (= Lange) 

Or in English: 

Snai IShel lX_L(b, n, t, tp) := SchneckenTorusX_L(b, M(b, n) , t, tp) 

Snai IShel lX_U(b, n, 0, tp) := SchneckenTorusX_U(b, M(b, n) , 0, tp) 

Snai IShel lXjAf(b, n, 0, tp) := SchneckenTorusX_W(b, M(b, n), 0, tp) 

Piotr provided a distance function which makes possible to plot three snail shells in one 3D plot: 


ABSTAND(tx, r) := r.[0, COS(o(), SIN(a)] 
DIST(ot, r) := r-[0, COS(a), SIN(ot)] 
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SchneckenHaLisXJ_(2, 20, t, s) + ABSTAND 

SchneckenHaiJsX_U(6, 20, t, s) + ABSTAND 

SchneckenHausX_W(e , 20, t, s) + ABSTAND 
Using the English function names: 

Snai lShel 1XJ_(2 , 20, t, s) + DIST 

SnailShell)UUC6, 20, t, s) + DIST 


■TT, 1 


TT, 1 

6 


+ ■ TT , 1 

. 2 


• n, 1 


■ TT , 1 


Snail She ll>UWCe, 20, t, s) + DIST 


+ 1 7T , 1 

. 2 



Very strange: One can plot the snails and enjoy their beautiful shapes. After saving the file 
as a dfw-file (as one is accustomed to) and then reloading plotting again does not work. One 
can find a strange message in the message bar (created by WRITE in function h(b,m,n). If 
saved as mth-file you will face no problems at all. Josef. (See also DNL#84 and the following 
comments.) 
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Comments and email-dialogue with Piotr Trebisz: 

My questions sent to Piotr: 

In the first figure I plot the spiral together with its respective anti-spiral and the corresponding 
tangents in their - trivial - intersection point. 


SPIRAL (2, 0.4, t) 

ANTI_SPI RAL (BASE (2 , 0.4), 0.4, t) 

[1, 0.4, 0] 

[1, 0.4, 0] + t ■ SPI RAL_TANGENT (2 , 0.4, 1) 

[1, 0.4, 0] + t ■ ANT I _S PI RAL_T ANGENT (BASE (2 , 0.4), 0.4, 1) 




The next figure should show the “surface ellipse”; but it is not an ellipse at all. 
ELLIPSES, 0,4, RADIUSJ<C2, 0.4), RADIUS_Y(2, 0.4), c (?) 



p 
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Piotr’s answers: 

To be honest, I am not very happy with your second figure. Graphs of both axes of the curvilinear 
coordinate system do not give a meaningful representation and impression of the whole system. I 
would prefer plotting a plane in this system. The parameter lines give a much better insight of the sys- 
tem. 

Example: b = 2 and m = 0.4 

X_ACHSE(2 , 0.4, t) + s ■ YJ\CHSE(2 , 0.4, t) 



The surface ellipse is an “ellipse” in the curvilinear coordinate system, not in the Cartesian one! 
This system is curved and has a torsion. This is the reason for the ellipse to appear as a twisted oval as 
shown in your third figure. 


My reply (and more questions): 

You are right with your representation of the coordinate system. I find the plot of the “twisted” 
ellipse very interesting and I will put it together with the “curved plane” and your comments 
into the next DNL. 

I have two more questions: 

(1) I wonder if it is possible to produce a longitudinal intersection of the shells. Is this 
possible using the curvilinear coordinate system? 

A nice view into the interior of the shell is given by parameter domains 0 < <f> < n (n/2) 
and 0 < t < 2 for the length-parameterized shells. I attach two plots. 

(2) I come back to a question posed in an earlier DNL, too. Why do you not use 
NSOLUTIONS for finding the zero of f(b,m)' ? Then you could do without the 
ITERATES- construct. 
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The interior of a shell = s): 



0 < t < 2, 0 <s<ji 



0<r< l,0<5<7r/2 


Piotr ’s next answers: 


Oh yes, it is possible to produce a longitudinal intersection. Solve the z-coordinate for 6 or t. This so- 
lution is say thetajcy. Then substitute theta_xy for 6 or t in the x- and y-coordinate. What remains is 
the intersection curve of the shell with the xy-plane. It is interesting that there are infinite many solu- 
tions of form thetajcy + A 2 ( k integer). Taking some As will result in a longitudinal intersection. 

I followed Piotr’s instructions: 


xy_t := (SOLUTIONS ((SchneckenHaus>LW_JB(e , 10, t, s)) =0, t)) 

3 

1 



r 

r k 

intersect := 

SUBST 

c 

SchneckenHausXJflLJB(e , 10, t, s) , t, xy_t + 

l 2 JJ 


[ 1 , 2 ] 

VECTOR (APPEND( -intersect, [0]), k, -2, 6) 
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The plot above is a top view (see the axes) and these are the front and side view. It is really 
the intersection with the xy-plane. 

z 


y 


X 


These plots and the next ones as well are produced with ...JB-function. I will explain later. 
See first some parameter curves of the (b = 2, m = 10) shell. 


VECTOR (SchneckenHausXJflLJB (2, 


VECTOR 


SchneckenHausXJflLJ B (2 , 


10, k, s), k, 1, 2, 0.05) 

77 ^ 

10, t, S), S, -77, 77, 

s ^ 



Piotr: 

What concerns function M(b,n) there are two reasons why I did not use NSOLUTIONS : 

1. I like to catch mathematical problems completely and to solve them by myself. Moreover my 
function is much faster because I defined the first derivative explicitly in a very compact form. 
And additionally there is a pretty annotation of the calculation progress (WRITE-command). 

2. As you found out the initial value 0 for m is a problem, because then thr derivative is 0, too 
which leads to a division by 0. My opinion is that fixing the lower bound of the solution interval 
for NSOLUTIONS as a small positive number is not a good idea because for small 6-values m tends 
to 0, i.e. for a too big lower bound NSOLUTIONS would not find a solution for very small 6-values. 
My short program uses the fact that j{b, m) starts increasing and then switches to decrease. I start 
with initial value 1 and increase until j{b, m) <0, then I turn and approach the zero from the other 
side applying the Newton-algorithm. Analysing f{b, m) shows that the zero of m < 4, i.e. finding 
the appropriate initial value does not need long time. 
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I accept Piotr’s explanations, see is my (Josefs) version using NSOLUTIONS. 


rn_(b, n) :z APPROX( (NSOLUTIONS (f (b , m) = 0, m, 0.1, «0) , n) 

1 

SchneckenHaus)CL_DB(b , n, t, cp) := SchneckeniTorLisX_L(b , m_(b, n), t, cp) 

SchneckenHausX_U_DB(b, n, 8, cp) := SchneckenTorusXUJ(b , m_(b, n), 8, cp) 

SchneckenHausX_W_JB(b , n, 8, cp) := SchneckenTorLisX_W(b , m_(b, n), 8, cp) 


The advantage of these JB-functions is that one does not face the problems when saving the 
file as dfw-file and then reloading it - as described on page 10. 

Final comment: 

Following Piotr’s DERIVE code you can see unusual variable names as: Ri, R 2 and m 0 . 

I wanted to find out how to create these variable names - it is not a subscript produced by 
R SUB 1 etc. 

Please follow my file CODE.dfw: 

#1: [Ri :=, R 2 :z > m o :z l 

#2: NAME_T0_C0DES(R 1 ) = [82, 127, 50, 48, 56, 49] 

#3: NAME_T0_C0DES(R 2 ) = [82, 127, 50, 48, 56, 50] 

I copied the variables in #1 from Piotr’s file and tried to find their “codes”. 

Ri and R 2 must be written under quotes as argument of NAME_TO_CODES. 

#4: C0DES_T0_NAME ( [ 82 , 127, 50, 48, 56, 49]) = R ± 

Ri appears as string (i.e. you can see “Ri” when copied into the edit line. 82 is the code of 
“R”. Then the remaining code numbers should be responsible for the subscript ‘V’? 

#5: N AME_T0_C0 DE S ( R ) = [82] 

#6: CODES_TQ_NAME ( [ 127 , 50, 48, 56, 49]) = ± 

#7: C0DES_T0_NAME ( [ 127 , 50, 48, 56, 50]) = 2 

The codes for "A" to "Z" are 65 to 90; the codes for "a" to "z" are 97 to 122. I try to produce 
variable names A^ a-i and z 40 . 

#8: CODES_TO_NAME([65 , 127, 50, 48, 56, 49]) =A ± 

#9: CODES_TG_NAME([65 + 32, 127, 50, 48, 56, 49]) = a ± 


#10: C0DES_T0_NAME ( [ 122 , 127, 50, 48, 56, 52, 127, 50, 48, 56, 48]) = z 40 
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Hi, it works. We have to remove the quotes and then we can use the newly created variable 


names: 


#11: 

[A x :z 10, ai :z 20] 

#12: 

Ai ■ a^ — 200 

#13: 

a :z 30 
1 

#14: 

A ± ■ ai ■ a — 200 ■ a 
1 1 


Please note the difference between a x and a 1 which is the result of a SUB 1 . 


I can not resist presenting two “non mathematical” snails: 




These are two snails which live on trees in order to escape their enemies. I found them in the 
woods of Madagascar, Josef. 
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Brussels Gate in Dendermonde, Belgium 

Dr.-lng. E. Lantschoot, Lantschoot_Weisel@web.de 
(1) The Program "gateQ", meaning the BRUSSELS GATE 

Since the end of the Middle Ages, the small Flemish city of Dendermonde has been a military 
stronghold. Situated at almost equal distance from Ghent, Antwerp and Brussels, it lies at the 
confluence of two rivers. Its system of defence was as simple as effective: by a controlled 
flooding of the entire surrounding area, invaders were kept away. Even the French king Louis 
XIV, in his endeavour of fortifying the north most boundaries of his territory, did not prevail. 
"Que n'ai-je une armee de canards...", "If I only had an army of ducks..." 

A monumental relict of those times is the Brussels gate. It was built in 1 822 under Dutch rule, 
and it featured a drawbridge which, when drawn, isolated the city. Although the drawbridge is 
now replaced by a fixed section and although the ropes disappeared, the counterweights (in 
the form of drums), the rails on which they rolled and the winches are still there. 



The item of interest to us is the profile of the rails. As it so appears, its shape was dictated by 
the intention of keeping the bridge and the counterweights in equilibrium. As we will see in 
the course of our discussion this observation is true, but only to a degree. 

From a mathematical point of view it is easier to first consider the theoretical trajectory of M, 
a point representing the center of the axis of the counterweight. This theoretical trajectory is 
found by applying the Runge-Kutta integration procedure rk23 (built in in TI-NspireCAS), and 
we see, - to our disappointment? - that it is not a quadratic curve, as we may have assumed. 
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Without being aware of Runge-Kutta, the Dutch military engineers deliberately replaced the 
theoretical trajectory by a theoretical quadratic curve as a surrogate. In doing this, a disequi- 
librium of forces was created in almost every point of the new trajectory. If this disequilibrium 
was positive for every point, the bridge would tend to lower itself until the support of the op- 
posed bridgehead was reached. To draw the bridge, a winch would have to be installed and 
it would have to develop a force at least equal to the disequilibrium force. 

The goals and the problems of a good design can be summarized as follows: 

1) To compute the theoretical equilibrium trajectory of M, with the weight of the drum 
as a parameter. 

2) To compute the surrogate quadratic curve - we call it the design-curve - with two 
fixed points and a third one, which influences its shape. 

3) To compute the disequilibrium force, occasionally rearranging 1) and 2) so that this 
force is positive, relatively low and relatively constant. 

4) To compute the total work to be supplied by the winch, and to find out whether the 
work can be done by one or two humans at each winch. 

5) With the data at our disposal to compute the profile of the wrought iron rails. We 

will see that it is also a quadratic, and since a quadratic can be designed by mechanical 
means, a wooden template could be constructed to serve as an aid in the forgoing 
process of the rails. 

In our programs, we are using row vectors because they fit better in a line of type. If one is 
more familiar with column vectors, it will help to consult the rules of transposition. 

A symbol denoting an initial status ends on 0 (zero). A force either begins or ends with an f. 
An arm, as used in crossp(force, arm) ends on an a, e.g. cga, center of gravity arm. 

Most data in our programs, such as the structural data length, height, weight etc. cannot be 
modified. Data by which the design can be altered, such as cga or adj, the adjusting factor of 
the weight of the drum, can be modified by the operator, but only in section gatel . 


We are using three important rules of statics: 

1) A force in a rope assumes the direction of the rope. If the direction of the rope is 
changed, e.g. by a deflection sheave, the direction of the force also changes, but its 
scalar value remains. 

2) A point on a curve is at rest, is in equilibrium, if the resultant of all forces acting on that 
point is perpendicular to the tangent of the curve at that point. To say it otherwise, the 
point is at rest if the resultant takes the direction of the normal. If the resultant has a di- 
rection deviating from the normal, then the projection of the resultant on the tangent is 
the force which will move the point. 

3) A rope is not a solid element, such as the girders of a bridge. It is treated as such as 
long as it is under stress, but when the rope is relaxed, it transmits no forces anymore. 
That is the reason why we required that the disequilibrium force should be positive 
(= traction, stress) for all states of the drum. 
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SECTIONS GATE1 AND GATE2 

The purpose of section gatel is to present, on the console, all actual data and to enable the 
operator, if he wishes so, to modify, either the arm of the center of gravity cga, or to adjust by 
adj the assumed weight of the drum. The assumed weight of the drum wdO is the weight 
which keeps the flat bridge in equilibrium. Why it can or should be adjusted is explained im- 
mediately. If the choices are made, the program proceeds automatically to section 2. 

First tool: the method of Runge-Kutta 

The initial position of the drum is precisely known: the distance between its center M and the 
deflection sheave U (= Umlenkrolle) (here reduced to a point) equals frr (free rope) and it 
hangs on the wall of the gate (distance M to wall: rd, radius of drum). The final position to the 
contrary is only defined by the fact that MU now equals lb*V2 + frr (rope at the bridge side 
now completely drawn to the drum side plus the length of the free rope), but it is not known at 
what final height above the street level. We expected a height of rd above null, but that is 
only wishful thinking; the actual height computed by Runge-Kutta for each of the 20 test 
points (including the final point M(20)) is dictated by the equilibrium of forces. By adjusting 
(with adj) the weight of the drum, we can, in the next run, either lower or lift the profile of 
Runge-Kutta. The range of adj is approximately between 0.95 and 1 .3. 

Second tool: the design curve 

If instead of implementing the Runge-Kutta curve we decide to replace it by a quadratic curve 
of our choice (thereby taking the risk of introducing disequilibria of forces), we can specify 
and we will specify the coordinates of M(1) and M(20). 

However, it takes three points to design a quadratic curve; we thus add a third point P3 p (m) 
to the right of M(20) and q (m) above it. In the Graphs-application p and q can be defined by 
shifting the sliders, and the shape of the design curve immediately adapts itself to the choice 
which was made. Caution however! The purpose of the design is not to closely adapt the 
shape of the design curve to that of Runge-Kutta, but rather to obtain a positive disequilib- 
rium force (see below). 

Program gatel (). (I don’t reprint the program code. What follows is the program output of 
program gatel, Josef) 

fOK/fll 
Locked data: 

Length of the bridge lb = 3.4 
Weight of the half bridge wb = 12000. 

Length of free rope frr =0.9 
Radius of drum rd = 0.25 

Data for beginner are: 

center of gravity arm cga = [-1.7 -0. 15] 

Adjusting factor adj = 1. 

Actual data are: 

center of gravity arm cga = [-1.9 -0.4] 

Adjusting factor adj = 1. 

new cga? (enter & OK/simply OK) [[-1 .9,-0. 4]] 

new adj? (enter & OK/simply OK) 1 . 

Unadjusted weight of drum wdO = 9483.55 
sol = 

[ 0.25 0.487395 0.724791 0.962186 1.19958 1.43698 1.67437 1.91177 2.14916 2.38656 2.62395 
[2.53542 1.81528 1.45212 1.17377 0.942654 0.744341 0.571291 0.418849 0.283963 0.164448 0.058665 H 
new design-curve? (y,n) n 









Program gate2() 


Define gate20= 

Prgm 

© Computes the true trajectory of M 

© Designs a quadratic curve, surrogate for the true trajectory 
Local autocodes JirJ\rv t rv t tv 
2 -» lastpr 
Lock cga } adj 

Unlock wdO: wd0\=force(wb } cga } 0 ): Lock wdO 
Disp "Unadjusted weight of drum wdO = " } wdO 
Unlock xb f xe,yb,ye 

xb\=rd © Initial location ofM } center of the drum axis 

xe\=j[j2 • lb+frr)^ ~[lb~rd)^ © Final location of M 

yb:=lb-jfrr 2 -rd 2 © Initial location of M 
ye\-rd © Final location of M 
Lock xbfXe^yb^ye 

© Computes, with Runge-Kutta, the true trajectory of M 
DelVar x } y 

fir\=force[yjb y cga>nQFrni\o ft]-[x y\)~fi'r) © Force in rope 
© fir as a vector firv, on the drum side ( direction has changed) 
firv:=fr- unitv([o lb]-[x y]) 
rv\=frv±adj' 0 -Wo] © Resulting vector 

tv\=rv 0 1 © Turning rv 90° to the left, has direction of the tangent line 

|-1 

yac :=—. © Inclination of the tangent line 

Mill] 

sol\=rk23\yac,x,y, {xb,xe \yb- 

\ 19 j 

Disp "sol = ",sol 
xlst : =mat ► 1 ist(so /_ 1 _ ) 
ylst : =mat ► 1 ist(so / 2 ) 

RequestStr "new design -curve? (y } n y\ant 
If ant="y" Then: 

© Design -curve 
DelVar p t q 

coJes:=rratMist(zeras({ a ■ xb^+b' xb+c-yb } a' xe^+b' xe+c-ye } a ■ [xe+p)^+b' (xe+p)+i 
des: =polyEval(c odes,x) 

Endlf 

EndPrem 


Program gate2 needs function force (see Appendix 1). 
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SECTION GATE3 

This section collects all the choices made in Section 2, and computes all items which are 
needed for visualizing the design in its present state. Section 3 has computed the vectors 
firv[n], force in the rope. These vectors have the direction from M to U, and that is consid- 
ered to be a positive direction, since it induces stress in the rope. If M(n) is not in equilibrium 
- and in most instances it is not - we must ask ourselves which mechanisms will cause M(n) 
to move. According to rule 2 of the statics, the moving force on M(n) is the projection of the 
resultant on the tangent. In a graphic construction, we depart from the tip of the resultant, 
take the direction of the normal, cross the tangent, proceed and cross the line MU, thus de- 
fining a point V. The section MV defines direction and strength of a vector, which will move 
M(n) when the system is left to its destiny. This vector is the disequilibrium force; if it has the 
positive direction (from M to V) it will lift the drum and lower the bridge. Therefore, section 3 
will warn the operator when the list of these vectors features a non-positive entry. Then, the 
design procedure of section 2 must be carried out again, but the following program permits 
us to visualize our mistakes first. 

(I don’t present program gate3(). I show the program’s output.) 

gales () r 

List of disequil. forces, diseqflst = 

{ 1950.82,2243.91,2234. 31, 21 14. 5, 1948. 46, 1763. 82, 1574. 86, 1389. 98, 1214. 63, 1052. 6, 906. 647, 778. 854, 670. j 
All the data needed for controlling the solution are now 
in storage. Please proceed with copt(n) and n of your choice 

Fertig 


THE PROGRAM COPT 

If one of the 20 test points is selected, copt(n) shows the diamond of forces in that point: 
force in the rope firv[n], weight of the drum adj*wd0, resultant r and the disequilibrium force 
diseqflst[n] as a part of force in the rope. In order to complete the picture, one can add the 
functions tali[n] (tangent in M(n)) or noli[n] (normal line) or even noli for the first 16 normal 
lines. The iron profile (computed in gate4) can be added later by entering the sets xilst and 
yilst. 


cop/(ll) 


Fertig 

diaf 

2.49378 

0.981556 



2.62395 

0.855318 



2.62395 

-1.04139 



1.86872 

-0.308971 



2.62395 

0.855318 



1.86872 

1.58774 



1.86872 

-0.308971 


xdiaflst 



{ 2. 49378,2. 62395,2. 62395, 1 . 86872,2. 62395, 1 . 

86872,1.86^ 








P 22 Erik van Lantschoot: Brussels Gate in Dendermode D-N-L#89 



SECTION GATE4 

This section provides: 

1 ) general data, such as potential energy when bridge high or low, and their difference; 

2) human data 

the above difference must be supplied by human work, and the human work is com- 
puted as the integral of diseqflst over the lenght of the rope drawn. Since diseqflst 
mostly is a list of steadily decreasing forces, it makes sense to supply these forces by a 
winch with a tapered drum (Spiral of Archimedes), on which an almost constant mo- 
ment cm can be applied. Finally, it is checked, what (human) force on the crank handle 
is needed. 
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3) the profile of the wrought iron rail. Since gate3 already provided the equations of all 
normal lines, this information can be used to compute a set of points l(n) which are rd 
below M(n); If a cubic regression function for the l(n) points is computed, it so turns out 
that the coefficient of x 3 is negligibly small. Consequently, the profile of the iron rail also 
is a quadratic curve. 


Output of program gate4() and respective plots: 


gate4[) 



pot. energy when bridge hi = 27600. J or Nm or Ws 

pot. energy when bridge lo = 21673.9 J or Nm or Ws 

pot. energy difference = 5926.11 J or Nm orWs 

For visualizing the human work, use the appL Data & Statistics 

Human work (compare wt. pot. en. diff.) = 5881.57 J or Nm orWs 

If to be supplied in 5 min, then power = 19.6052 W 

Force at the crank handle = 167.157 N 


For visualizing the wrought iron rail, use the appl. 
Graphs and enter the lists xilst and yilst 

Per tig 

I of cm 

1459.95-98.4747- x 

i 0 

2/99 
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(2) DERIVATION of the Function “Force” 

The statement of the problem is simple. A bridge 
with a hub at (0,0) and an end at (-lb,0) should be 
hauled up by means of a rope which runs over a 
deflection sheave at (0,1b). The weight of the 
bridge is wt N, and the arm of the gravity center 
is a vector [r,s] (both elements are negative). 
The problem is to express the force in the rope 
as a function of the length of rope drawn. It’s a 
horrendous exercise in trigonometric fitness! 

Note - Because in this document "force" is 
treated as a standalone, we converted it to a 
program and commanded "disp" instead of "re- 
turn". 


The whole computation relies on equating the 
moment of the weight and the moment of the 
force in the rope. When moments are consid- 
ered, one is allowed to lenghten or to shorten an 
arm, provided the force acting on that arm is 
reduced or increased correspondingly. Here we 
lengthen the arm of the gravity center until its first 
component equals lb. The moment of the weight 
now is: 

mw = (w*r/lb)*norm([r*lb/r s*lb/r])*sin(...) 

... being the angle between the weight vector and 
its arm. 



Figure 1 


Figure 2 
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Figure 3 shows an equivalent situation: the arm 
of the gravity center has been decomposed in 
two vectors OH and OJ and the weight vector 
now acts on the points H and J of the bridge. The 
first moment is: 

ml = (wt*r/lb) *lb*sin(9O°-0) 

After simplification: 


We now pay attention to the vector OJ and 
lengthen it until its length equals lb. The moment 
then becomes: 

m2 = (wt*r/lb)*(s/r)*(s*lb/r)*(r/s)*sin(18O°-0) 
After simplification: 
m2 = wt*s*sin(0). 


ml = wt*r*cos(0) 

and the length of the other arm, OJ, equals 

s*lb/r. 



We now consider the scalar variable fir, force in 
the rope. We subdivide the segment H-SH in two 
equal parts H-M and M-SH and denote the angle 
H-O-M by p. 

mfir = fir*lb*sin(90°-p) = fir*lb*cos(p) 

and the general equation of the system obviously 
is 

mfir+m1+m2 = 0 

Our attention now goes to the line segment 
H-SH. If the amount (length) of rope with- 
drawn is denoted by xt , the length of H-SH 
is 



Figure 4 



V2*lb-xt and for the rectangular triangle H-O-M we find that the quantity w 

w = (V2*lb-xt)/(2*lb) = sin((90°-6)/2) = sin(p). 

We now must go back to the trigonometric functions in ml, m2 and mfir. 

In ml, we have sin(9O°-0) = sin((9O°-0)/2)= sin(2*p) = 2*sin(p)*cos(p) = 2*wW(1-w A 2). 

In m2, we have sin(0) = cos(2*p) = 1-2* (sin(p) A 2) = 1-2*w A 2. 

In mfir we have cos(p) = V(1-w A 2). 

The rest is a little bit of calculus. 

The advantage of this derivation is to reveal a detail which goes unnoticed when using 
crossp's: if the center of gravity is on the surface of the bridge, the force fir decreases line- 
arly with xt . 
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(3) The program mechtor() 

The gate to Mechelen (background picture on page 3 of the document) is an exact replica of 
the gate to Brussels. The approach we take her is that of an ordinary but mathematically in- 
terested tourist who discovers this contraption, takes a few pictures, notes a few dimensions, 
and then tries to find the mathematical equation of the iron rail by means of the here pre- 
sented program mechtor(). 

When looking at page 3, we see a background picture of the Mechelen Gate and a function 
graph as an overlay; the graph is as left by the last user of the program. 



Therefore, we must first explain how the graph was originally designed and what it means. At 
the very first use of the program, a list of x-values xlst:=seq(n-1,n,1,20), the list of numbers 
0 to 19, was created. This list was locked. Then a second list ylst := newlist(20), not locked, 
was defined. Using the Data & Statistics Application the 20 points (xlst[n], ylst[n]) were dis- 
played as a horizontal row of red points. A movable line m2(x):=a + b*x was added; this line 
the footprint of the rail is supposed to be the vertical projection of the iron rail on the ground. 

These raw data are now to be moved to their proper position in the picture (see the respec- 
tive instructions in the TI-Nspire_SS_Guide, English, pp. 334 and German pp. 374). First of 
all the origin must be moved until the y-axis passes the highest point of the rail. Then the x- 
axis must be stretched or condensed until xlst[20] lies exactly under the point where the rail 
meets the street. 

The red points El to E20 with En = (n-1,0) are now moved to their positions on the rail by 
grasping them one after the other with the cursor, and moving them up or down (moving right 
or left is excluded because xlst has been locked). 
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Note: a point which is grasped loses its color. When the point is on its proper place then click 
on it again, and it recovers its color. The reason is, that the movement of the cursor is ap- 
plied to all non colored points . Finally the footprint is grasped and dragged to its proper po- 
sition. Then mechtor() was continued and the 4 th degree regression curve fotreg of the 
points E1-E20 was calculated. The white curve displayed on page 3 is fotreg (if its equation 
is not shown then simply click on the curve). This curve, however is not the solution of our 
problem, it is only the key for it. fotreg indeed is only the regression curve approximation 
of the photographed iron rail, and must be corrected for effects of perspective and scaling. 

What can we do now, which has not been done already by the last user? 

On page 3 - and without opening the program - we can move one or the other point on the 
rail, thereby observing how fotreg changes. We can also move the footprint. We call these 
adjustments “corrections”. 

If we would like to know the equation of the rail in metric units, with the origin of meas- 
urement (0,0) at the intersection of the footprint with the wall then we have to run the pro- 
gram by entering mechtor() on the calculator page. 

Immediately appears the question “new data or corrections? (nd/c)”. If the answer is “c”, then 
the computation is performed with the most recent data given on page 3; if the answer is 
“nd”, then points El to E20 and the footprint as well are set back to some initial format and 
the whole process bringing everything to the proper place must be done again. 

Now we have to explain how the program converts fotreg to the final answer rail_c which is 
the equation of the rail in reality i.e. in metric units. The conversion needs four steps: 

a) We need to know the real measures. By inspection on the site we find El =(0,2.25) and 
E20=(4.82,0). 

b) The y-values of E1-E20 should not refer to the x-axis of the diagram but to footprint. 
The result of this correction is rail_a. 

c) rail_a will now be converted to rail_b in such a way that rail_b=0 for xlst[20]=4.82. For 
reaching this goal the list of coefficients cof of rail_a are multiplied by list caf. 

d) The ordinates of rail_b will now be multiplied by such a factor that ylst[1]=2.25. 

The final result - of all this squeezing and stretching) is rail_c and it is left to the reader 
whether the coefficients of x 4 and x 3 are small enough to be neglected, thus supporting the 
statement that the equation of the rail is a parabola. 


mechtor [ ) 

new dat 

3, or corrections only? (nd,c) c 



fotreg = 

0.000002' x 4 -0. 000078- 

x 3 +0. 002697- 

x 2 -0. 084425- x+0. 926316 


rail_a = 

0.000002- x 4 -0. 000078- 

x 3 +0. 002697- 

x 2 -0. 077395- x+0. 772516 


rail_b = 

0.000489' x 4 -0. 004816- 

x 3 +0. 042253- 

x 2 -0. 306355- x+0. 772516 


rail_c = 

0.001423' x 4 -0. 014028- 

x 3 +0. 123066- 

x 2 -0. 892277' x+2.25 


Fertig 

rai/|.x-0 




2.25 

roi/|x=4. 8 

2 


0 

.005709 
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Comments of the editor: 

I do not reprint the program in order to save space. You can find it among the downloadable files together with 
the picture of the Gate. Thanks to Erik van Lantschoot for providing this. 

As the DNL-readers will already know I am an enthusiastic user of background pictures for several reasons. I 
like this example very much because it is much more than only modelling the picture of the rail; it does an im- 
portant step further: converting the picture to reality. I studied Geometry and I remember “Photogrammetry” 
where we had to reconstruct buildings from photographs considering perspective distortions. This gate is a great 
example what can be done with students (neglecting the perspective because being so close to the object). 

I had some very fruitful email exchanges with Erik van Lantschoot and in the last one I asked him why not im- 
mediately performing a quadratic regression? (You are invited to do this and then repeat the conversion to real- 
ity.) I liked his answer: “I assume that a tourist who has such a powerful tool like TI-Nspire at his disposal will 
choose the “best” available regression.” 

My experience as teacher is that most of the students think: “All curves are parabolas.” 

I found a pretty picture of Dendermonde from earlier times. Both gates mentioned in the contribution can be 
found there. I sent this picture to EvL. He wrote back: 


I lived in the centre of the town where the bridge crosses the Dender River until 1945. This was the 
year of my end examination at the gymnasium, then I had to leave for starting my studies at university. 



41. 'Porta XLechCiniensis, 42. Porta PruxeCCensis 


You may now find out EvL’s age of life - and he is doing mathematics with modern technologies. This is really 
great, congratulations. 
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Structures and Rules of the Prime Numbers 

(Observations and Conclusions) 

Dietmar Oertel, Germany 


Preliminary Remarks 

This paper deals with the phenomenon of the prime numbers. 

My observations with respect to the prime numbers are 

• structures in the set of primes can be observed if they are not considered all together in a formula 
but one after the other 

• a regularity in form of a delta-sequence with a certain sequence width appears which depends on 
the number of the considered prime numbers 

• this sequence width related to the given prime numbers and the number of the remaining indivisi- 
ble numbers can easily be found 

• the density of the prime numbers is decreasing 

• referring to the number of primes considered equations can be “constructed” which - based on 
combinatorics - can determine the n+l th prime number if n primes are known 

• understanding the delta sequence and its width it becomes clear that there is no regularity in the 
sequence of the prime numbers when viewing all primes at a whole 

• my results match the Euler Equation (I learned this in a TV broadcast from 5 September 2011 
‘hitec’ where this equation and the Zeta function were presented). 

http://www.3sat.de/page/?source=/hitec/156496/index.html 


My summary of these observations: 

It is not difficult to recognize the structure of the prime numbers considering the following aspects: 

* prime numbers shall one after the other step by step be considered in an expandable formula 

* the principle of combinatorics based on Pascal’s Triangle is applied 

* only the integer parts (FLOOR- function) of possible fractions appearing in the formulae must be 
used 


My final result is the following formula (MathCad and DERIVE): 

v p) ! 




"Max 


langsj 

SPrmFunPNI 

■TV n Max) : 

= I 

ti = 0 

floor(n) + 

i 

k = 


nMax 

f 

DIM(vPrm) 

k d: 

SPrmFunPN(vPrm, 

nMax) := I 
n=0 

FLOOR(n) 

+ I 

k=l 

(-D ■ 


k! ■ ( langej J - kj 1 


i 

j = 0 


I 


floor 


I K PrinTl I v , kl I 

i 


FLOOR 


(KmbNumPrdRev(vPrm , k)) 


J JJ 


(The principle of combinatorics is represented in KPrmll, which will be explained later.) 
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We will see that the next prime number is found when SPrmFunPN(v Prm ,k) = k + 1 or, defined as a 
function: 

PrnNul(v p ,k) = SPrmFunPN(v p ,k) - 1 - k. 


Example: 

v Prm 0 := 2, v Prm := 3, v Prm := 5; v Prm r =(2 3 5) the next prime number follows if the following condi- 
tions are fulfilled: 

PrnNul(v Prim , k) = 0 and k > 0 leads iok=l and v Prm3 := 7 . 


I used two Mathematics programs for my investigations, MathCad 14 and DERIVE 6.10. 

DERIVE 6.10: Advantages; very fast and adjustable accuracy 

comfortable building functions 
Disadvantage: DERIVE has been taken off the market 

MathCad 14: Advantages: comfortable presentation of functions 

Disadvantage: restricted accuracy 


Regularities and Delta Sequence 


As I noted above it is not possible to recognize any regularity in the distances between the prime num- 
bers if their sequence is investigated as a whole. But it is possible performing a step by step approach. 
Some graphics shall demonstrate this in the following. 

Take the natural numbers (without 0) then the numbers indivisible by 2 are the odd numbers and their 
differences form the “Delta Sequence” {2, 2, 2, 2, 2, ...}. 

Remove the numbers divisible by 2 and 3 (base set {2,3}) from the natural numbers then the remain- 
ing numbers are {1,5,7,11,13,17,19,23,25, ...} which results in the delta sequence {4, 2, 4, 2, 4, 2,...}. 

Base set {2, 3, 5} gives a delta sequence {6, 4, 2, 4, 2, 4, 6, 2, 6, 4, 2, 4, 2, 4, 6, 2, 6, 4, 2, 4, 2, 4, 6, 2, ...} 

You will recognize the periodicity the sequences. 

Will we find this periodicity for bases sets {2, 3, 5, 7} and {2,3,5,7,11} and {pi,p 2 ,p 3 , ...,pk}, too? with 
Pi the i-th prime number? 

Can we predict the period length of the delta sequence? 


See two graphic representations for base sets {2,3,5} and {2, 3, 5, 7}: 


ei 

19 

17 

15 

13 

11 

9 

7 

5 

3 

1 


I I 


I I I 


ro ro ro pj ro u to to to to 


ui tn lh c_n 



Width, Number and Sum 

Note (Josef): I wanted to illustrate the delta sequences using DERIVE and defined four func- 
tions: 

nd(n, k, p, st, i) := 

Prog 

st := VECTOR (k_, k_, 1, n) 

i := 1 
p := 2 

#1 : Loop 

If i > k 
RETURN st 

st := SELECT (M0D(v, p) t 0, v, st) 
p := NEXT_PRIME(p) 

i :+ 1 

n_d(n, k, p, st, i) := 

Prog 

st := VECTOR (k_, k_, 1, n) 

i := 1 
p 2 

#2: Loop 

If i > k 

RETURN VECTORClv, 1], v, st) 
st := SELECT (M0D(v, p) t 0, v, st) 
p := NEXT_PRIME(p) 

i :+ 1 

nd(n,k) returns the sequence of all integers 0 < x < n which are indivisible by the first k 
prime numbers. 

n_d(n,k) marks all these integers from nd(n,k) as a point (x,1) 

Examples: 

#17: nd(300, 4) 

#18: [1, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 

79, 83, 89, 97, 101, 103, 107, 109, 113, 121, 127, 131, 137, 139, 143, 
149, 151, 157, 163, 167, 169, 173, 179, 181, 187, 191, 193, 197, 199, 


209, 211, 221, 223, 227, 229, 233, 239, 241, 247, 251, 253, 257, 263, 
269, 271, 277, 281, 283, 289, 293, 299] 
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#21: n_d(300, 4) 


50 100 150 200 250 300 

(sequence of all integers < 300 which are indivisible by 2, 3, 5, and 7) 

#23: n_d(300, 5) 


50 100 150 200 250 300 

The periods are not really clear to recognize. Two more functions shall help: 

d_seq(n, k) := VECTOR((nd(n, k)) - (nd(n, k)) ,1,1, 

#3 : i + 1 i 

DIM(nd(n, k)) - 1) 

A_seq(n, k) V ECTOR (Ti , (n_d(n, k)) - (n_d(n, 

#4: L i + 1,1 

k» 1,1,1, DIM(n_d(n , k)) - 1) 

i,lJ 

d_seq(n,k) gives the delta sequence of nd(k): 


|d_seq(50, 2) 


[4, 2, 4, 2, 4, 2, 4, 2, 4, 2, 4, 2, 4, 2, 4, 2] 


|d_seq(100, 3)| 


[6, 4, 2, 4, 2, 4, 6, 2, 6, 4, 2, 4, 2, 4, 6, 2, 6, 4, 2, 4, 2, 4, 6, 

2 , 6 ] 


|d_seq(300, 4) 


[10, 

2, 

4, 

2, 

4, 

6, 

2, 

6, 

4, 

2, 

4, 

6, 

6, 

2, 

6, 

4, 

2, 

6, 

4, 

6, 

8, 

4, 

2, 

4, 

2, 

4, 

8, 

6, 

4, 

6, 

2, 

4, 

6, 

2, 

6, 

6, 

4, 

2, 

4, 

6, 

2, 

6, 

4, 

2, 


4, 2, 10, 2, 10, 2, 4, 2, 4, 6, 2, 6, 4, 2, 4, 6, 6, 2, 6, 4, 2, 6, 

4, 6] 

This is a little bit better but still not very illustrative. A_seq(n,k) plots the points of the se- 
quence one after the other, like (1,10), (2,2), (3,4), (4,2), ... for the 4-delta sequence from 
above. Let’s inspect the graphs for k = 2, 3, 4, and 5: 

As I noted above it is not possible to recognize any regularity in the distances between the 
prime numbers if their sequence is investigated as a whole. But it is possible performing a 
step by step approach. Some graphics shall demonstrate this in the following. 






A_seq(300, 3) 




Do you really see the periods? Compare with the nd-lists and the n_d-graphs from above! We 
will proceed for k = 5: 


#11: n_d(300, 5) 


-l 


50 


100 


150 


200 


#12: A_seq(1000 , 5) 



250 


We cannot discover any period. Maybe that we need a larger n? Take n = 2000: 
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Obviously there seems to exist a symmetry. 
k = 3, n = 300, period length = 8 

[ 6 , 4 , 2 , 4 , 2 , 4 , 6 , 2 , 6 , 4 , 2 , 4 , 2 , 4 , 6 , 2 , 6 , 4 , 2 , 4 , 2 , 4 , 6 , 2 , 6 ,...] 


k = 4, n = 500, period length = 48 

[ 10 , 2 , 4 , 2 , 4 , 6 , 2 , 6 , 4 , 2 , 4 , 6 , 6 , 2 , 6 , 4 , 2 , 6 , 4 , 6 , 8 , 4 , 2 , 4 , 2 , 4 , 8 , 6 , 4 , 6 , 2 , 4 , 6 , 2 , 6 , 6 , 4 , 2 , 4 , 

6 . 2 . 6 . 4 . 2 . 4 . 2 . 10 . 2 . 10 . 2 . 4 . 2 . 4 . 6 . 2 . 6 . 4 . 2 . 4 . 6 . 6 . 2 . 6 . 4 . 2 . 6 . 4 . 6 . 8 . 4 . 2 . 4 . 2 . 4 . 8 . 6 . 4 . 6 , 

2 . 4 . 6 . 2 . 6 . 6 . 4 . 2 . 4 . 6 . 2 . 6 . 4 . 2 . 4 . 2 . 10 . 2 . 10 . 2 . 4 . 2 . 4 . 6 . 2 . 6 . 4 . 2 . 4 . 6 . 6 . 2 . 6 . 4 . 2 . 6 , .. .] 


k = 5, n = 10 000, period length = ??? 

According to the graph from above we can assume that it is about 500. 


Back to Dietmar’s paper: 

I will define a width as follows: 

width (k) is the number of all integers containing the differences of one period, or said with other 
words, it is the first divisible number after the first set of indivisible numbers by the first k prime num- 
bers. This seems to be a bit complicated. I will make this clearer using a short table: 


k 

base set 

width 

number = 
period length 

indivisible integers 
(including 1) 

sum of the 
indivisibles 

1 

{2} 

2 

i 

{i} 

i 

2 

{2,3} 

6 

2 

{1,5} 

6 

3 

{2,3,5} 

30 

8 

{1,7,11,13,17,19,23,29} 

120 

4 

{2, 3, 5, 7} 

210 

48 

{1,11,13,17,19,23,29,31,37 

? 

41,43,47,53,59,61,67,71,73 

79,83,89,97,101,103,107, 

109,113,121,127,131,137, 

139,143,149,151,157,163, 

167,169,173,179,181,187, 

191,193,197,199,209} 

5040 

5 

{2,3,5,7,11} 

2310 

480 




General formulae are: 


k k 

width(&) = P[ p t , number (/() = ri(A-i) 

i=\ i = 1 

According to these conjectures (there are no proofs) the period length for k = 5 should be 480. The 
plot seems to confirm this. I will check it supported by DERIVE: 
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setl := (d_seq(3000, 5)) 

[1 480] 

set2 := (d_seq(6000, 5)) 

[481 960] 


setl - set2 


The result is a list of 480 zeros. So I should be right. 

What concerns the sum of one period it is remarkable that the sums of the first and the last, of the sec- 
ond and the last but one, of the third and the ... are constant and equal the width. 

e.g. for A: = 4: 1 +209= 11 + 199 = 13 + 197 = ... = 103 + 107 = 210 (24 times) 

This makes it easy to calculate the sum: 

j 1 k k 

sumper(£) = — • width(£) • number(^) = ~ Y\ Pi n (Pi ~ 1) • 

2 2 i=l i=l 

Let’s sum this up with DERIVE: 
width(k) := n(VECT0R(NTH_PRIME(i) l i, 1, k)) 
numb(k) := II(VECT0R(NTH_PRIME(i) - 1, i, 1, k)) 

VECTOR (width(k), k, 1, S) = [2, 6, 30, 210, 2310, 30030, 310310, 9699690] 

VECTOR (numb (-i), i, 1, S) = [1, 2, S, 4S, 4S0, 3760, 92160, 163SSS0] 


width(k) ■ numb(k) 



VECTOR (sumper(-i), i, 1, S) 

[1, 6, 120, 3040, 334400, S64S6400, 23324300S00, S043310S73600] 


I will define a kind of “prime number density” with respect to the number of the first k prime numbers: 

numb(£) 


my_pr_dcns(/v) = 


width(^) 


I know that not all numbers in numb(k) are prime numbers because only the divisibility by the first k 
primes is considered but however we can see that this “density” also decreases. 

my_pr_dens(x) 


VECT0R([x, my_pr_dens(x)] , x, 1, 100) 

■ &-.5 


0.4 


0.3 



10 20 30 40 50 60 + 70 SO 90 100 
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0.5 


2 4 6 S 10 12 14 16 IS 20 


The plot of my_pr_dens(x) appears as a piecewise defined function. 

This “density” must not be mixed up with the prime number density function tt(x) how it is 
well known in mathematics. tt(x) gives an approximation of the number of primes < x. 
tt(x) was - as far as I do know - defined by Gauss as follows: 


tt(x) = 


— = Li(x) (= Logarithmic Integral). 
J Inf 


2 


Using this counting method we can define the density as 


Li(x) 


Dietmar mentions the “Euler equation” which was discovered by Leonhard Euler: 

1 2 2 

Pi _TT 2 


This equation reads as lim FT — = — . 

“-*"1=1 Pi -1 6 


Again Dietmar Oertel: 

Similar to the numb(k) function from above I define a function numb_pl us(k) : 

numb_plus(k) := Tl(VECTOR(NTH_PRIME(i) + 1, i , 1, k)) 

VECTOR (numb_pl us (i) , i, 1, 8) 

[3, 12, 72, 576, 6912, 96768, 1741824, 34836480] 

But unlike to numb(k) and width(k) I cannot find any meaningful interpretation of these numbers 
except this one, that numb_pl us is necessary to connect my functions with the Euler equation: 

2 

width(k) 

my_euler_(k) 

numb(k)*numb_plus(k) 

my_euler_(100) = 1.644515221 

my_euler_(1000) = 1.644913174 


It needs a calculation time of 24 seconds for k = 1000. 
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my_euler(k) :r TT 


VECTOR 


2 

NTH_PRIME(i) 

2 

NTH_PRIME(i ) - 1 




k 


j j 


my_euler(10) = 1.63307049 
my_euler(100) = 1.644515221 
my_euler(1000) = 1.644913174 

This procedure needs only 16 seconds for k = 1000. 
Interestingly my third procedure is the fastest one: 


f_pr(k) := VECTOR (NTH_PRIME(-i), i , k) 

DIM (vprm) 

TT (vprm - 1) 
n=l n 

prm_"in1:_dens(vprm) := 

DIM (vprm) 

TT vprm 
n=l n 

DIM (vprm) 

FI (vprm + 1) 
n=l n 

prm_int_dens_p"lus(vprni) := 

DIM (vprm) 

TT vprm 
n=l n 


1 

prm_euler(vprm) := 

prm_"int_den3(vprm) ■ prm_int_dens_plus(vprm) 

prm_euler_f i rst(n) := prm_euler(f_pr(n)) 

prm_euler_f i rst(10) = 1.63307049 

prm_euler_f i rst(100) = 1.644515221 

prm_euler_f i rst(1000) = 1.644913174 

prm_euler_f "i rst (10000) = j] 

It needs only 8 seconds for k = 1 000, but 986 seconds for k= 10 000. 

2 
77 


. 644932811 


6 


= 1.644934066 
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Dietmar was inspired by my graphics (pages 33 & 34) that was eager to produce his own 
illustrations of the A-sequences periodicity. 


See two of his very pretty plots. 





(will be continued) 

Links: 

http://www.austromath.at/medienvielfalt/materialien/krvpto/lempfad/content/k prim.htm 

http://wwwmath.uni-muenster.de/u/deninger/about/seminare/vorlesungsskript-an- 
zahltheorie/primzahlen2 .pdf 

http://arxiv.org/ftp/arxiv/papers/0803/08Q3.042Q.pdf 

http://sciblogs.co.nz/guestwork/20 1 0/07/1 9/unexpected-prime-number-breakthrough/ 


p 40 


Derive- and CAS-TI-User Forum 


D-N-L#89 


Dear Josef, 


it's been such a long time ! ! ! 

I see you are well, always busy. Well done! I also try to keep fit, both physically and mentally. 

I remember from old times that you were very keen on 2D graphic representations in Derive. 

I’m also a fan of them, and I have recently come across Delaunay triangulations and Voronoi dia- 
grams. 

They are truly amazing and beautiful, and I'd like to have a Derive program to generate them, but 
I do not quite understand the algorithm. Can you help me? maybe you have developed a program 
yourself, I 

wouldn't be surprised if you did! ! 

thanks in advance, old friend 

Enric Puig, enricpuig026gmail . com 


DNL ; I must admit that I had never heard before about Delaunay (Voronoi sounds better for 
me but it is a very “weak” sound ...). I got some information from WikiPedia: 

A Delaunay triangulation for a set P of points in a plane is a triangulation DT(P) such that 
no point in P is inside the circumcircle of any triangle in DT(P). Delaunay triangulations 
maximize the minimum angle of all the angles of the triangles in the triangulation. 

The Delaunay triangulation of a discrete point set P in general position corresponds to the 
dual graph of the Voronoi tessellation for P. 




Dear Mr. Bohm, 

I read DNL#88 with much interest. I specially liked the contribution about T- and Z-scores. I’d like to 
add a note from practice. Statistical evaluations like these are made in medicine for osteoporosis pa- 
tients when measuring their density of bones. 

X-raying the lower lumbra vertebrae (untere Lendenwirbel) it is noticed how much of the given ray 
energy interfuses the observed bone. When much of the energy is absorbed then the bone is dense and 
stable, otherwise this is a sign for a weak bone with possible danger to be broken. 

Due to his age of life the bone density of everyone is decreasing and there is a graphic correlation 
between the age of a patient and his/her “age-related normal” bone density. Deviations are expressed 
by T- and Z-score. 
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At the moment I cannot write more about this but I will inform myself in the next future and then 
I will give the respective report. 

I have two more problems to discuss: 

1. ) In my materials from my time as student I found calculations for long division of polynomials. 
These were exercises in the 1 st semester. I still know how to perform a long division but I never 
heard for which purpose it is needed? In the meanwhile I found out that it can be used to reduce 
the order of a polynomial equation by dividing by a linear factor (x - solution) if a solution of the 
equations is known. This is clear for me. But in my old materials I had to divide a polynomial of 
order five by a cubic polynomial. I cannot realize any sense in this calculation! 

2. ) In my job I often had to deal with the problem finding an appropriate function describing 
physical measuring results. Since the EXCEL era this no longer a problem because after some tries 
a satisfying fit- function can be found. 

But there are physical processes which cannot be described by polynomial functions. In oil 
hydraulics one can measure spring and damping properties of plain bearings (= Gleitlager) and 
then present the results in a diagram. These diagrams are used for calculating the resonance behav- 
iour of fast running shafts (turbines or compressors). Because lack of a matching “polynomial ap- 
proximation” the curves are represented by many nodes for applying linear interpolation between 
them. This is not so bad but I’d prefer using an approximating function. 

In the book “ Bronstein-Semendjajew“ I found a type of curve which seem to fit - at least opti- 
cally 64 . It has the form: y = a + b/x + c/x 2 . 

Is it possible to find the coefficients a, b, and c supported by many nodes (xi/yi)? 


Many thanks in advance and best regards from sunny but cold Nuremberg, 
Klaus Korner, koerner . klaus6gmx . de 


DNL: Many thanks for your comments on T- and Z-scores. I am quite sure that Guido 
Herweyers will be very happy reading your kind words. 

What concerns your two questions: 

(i)i can imagine that students who find long division hard and boring might ask the teacher 
about its sense. What are possible answers? There are some applications of this algorithm, 
most of them beyond the students’ mathematics knowledge at the time when it is in their cur- 
riculum. 

• Partial fractions decomposition 

• Reduction of the order of a polynomial equation by a linear factor division 

• Investigation of the behaviour in infinity of rational functions 

• Calculating the proof sum in the CRC polynomial 

• One important step in the AES-algorithm (Advanced Encryption Standard) 

• GCD of two polynomials (Euclidean algorithm) 

• Function value, tangents & more for polynomial functions 111 

• possibly in the theory of Groebner Bases (I am not quite sure) 

111 Strickland-Constable, Charles, "A simple method for finding tangents to polynomial graphs", Mathematical 
Gazette 89, November 2005: 466-467. 
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I would like to give some examples: 

Investigate the behaviour in infinity for the given functions: 


... . x 4 - 3x 3 - 2x + 1 ... . x 4 -3x 3 -2x + 1 ... . x 4 -3x 3 -2x + l 

/!(■*) = t 3 , / 2 (x) = —z , / 3(x) = - 


2x -x 


2x -x 


2x-l 


We perform the division. The rational part tends to 0 for x ^ The integer part is the equa- 
tion of the asymptote (curve): 

4 3 


x - 3-x - 2-x + 1 


3 

2 «x - x 


5-x 


7 


2 

2- (2-x - 1) 


2 

2- (2-x - 1) 






We have a straight line, a parabola and finally a third order polynomial appearing as asymp- 
totes. In DERIVE we can receive the integer part and the remainder of the result immedi- 
ately. 

TI-NspireCAS offers polyQuotient and polyRemainder together with expand and propFrac. The 
V200 has only the last two functions available. 
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EXPAND 


4 3 

x + 3 -x + 3 -x + 6 


x - 5 


20-x 31 2 

+ + x + 3«x + 3 


2 2 
x - 3 x - 3 


4 3 2 2 

QUOTIENT (x + 3-x + 3-x + 6, x - 3) = x + 3-x + 3 

4 3 2 

REMAINDER (x + 3-x + 3-x + 6, x - 3) = 20-x + 31 


polyQuotient(x 4 -3' x^-2- x+1,2- x 2 -xj x 2 5-x 5 


2 4 8 


po ly Rema inder (x ^ - 3 ■ x ^ ! 

2- x+1,2- x 2 -xj i 21 x 

8 

expand | 

x 4 -3' x 3 -2' x+1 j 

-5 ^x 2 5x 1 5 

2-x 2 -x J 

8- (2- x-l) 2 4x8 

propFrac 

|x 4 -3' x" — 2 ■ x+1 

-(21-x-s) | x 2 5-x 5 

\ 2x 2 -x 

j 8- x (2' x-l) 2 4 8 


This is to tell about the integer part, but there is also a nice feature of the remainder. 



This seems to be the tangent, but in which point? What we do know is that dividing by a lin- 
ear factor with a zero gives the reduced polynomial and dividing by another linear factor 
(x - a) gives the function value f(a) as remainder. 

Both procedures can be substituted by the Horner algorithm - so it does not need long divi- 


sion. 
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Now the point is fixed! 

What about proceeding, i.e. let us divide by (x - a ) 3 and plot the remainder. What do you 
expect as outcome? 


3 


REMAINDER 


X 


L 2 


2 

2-x 

+ x - 3, 

5 


r 


x 


1 V 

2 J ^ 


2 

7 «x 5-x 47 


20 S 16 



It is the osculating quadratic function - a parabola. 



3 2 

x 2 «x 


TAYLOR 


L 2 


+ x - 3, x , , 0 

5 2 


203 


SO 


^3 2 

x 2 «x 


TAYLOR 


L 2 


+ x - 3, x, , 1 

5 2 


39- x - 121 


40 


^3 2 

x 2 «x 


TAYLOR 


L 2 


+ x - 3 , x , , 2 

5 2 


2S-x + 30 -x - 233 


SO 


4 3 2 3 2 

TAYLOR (3 -x - 14-x + x - 3 -x + 3 , x , 2 , 3) = 10-x - 71-x + 91-x - 43 


This does not need any comment. 
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Possible question for students: 

Do you have an idea how to extend Horner’s Rule in order to find the remainders from above 
without applying polynomial divisions? 


(2) To your second question: 

It is possible to find the regression line of the 
requested form. In DNL#63 I presented a gen- 
eralised regression procedure based on the 
normal matrix. 

DERIVE has this algorithm built in in its FIT- 
function: 

"1 -4 " 

2 3 



3 5 


data 


5 5 
26 

6 

5 


27 

7 

5 



f 

b c ' 


FIT 


x , a + + 

, data 



x 2 




X 

j 


In DNL#63 I showed a small program per- 
forming this regression: 



0.7264704722 


X 


10 . 1964S749 

+ 5.4459537S7 

2 

x 


And in DNL#89 I present the Nspire-form. 

I am using lists for x- and y-values because 
then plotting the scatter diagram is easier. 
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It might be a nice experience for student to perform the Euclidean algorithm for finding the 
GCD for polynomials. They can practise long division and then check their manipulations 
using the respective CAS functions. 

Let a = x 5 -2x 4 +2x 2 -13x+18 and b = x 3 -x 2 +8x-20. Find the GCD(a,jfc>). 

I present the manual procedure follwed first by the DERIVE treatment and then by the 
TI-NspireCAS screen shot. For obtaining the GCD only the remainders are important. But for 
controlling the long divisions the students would also need the quotients. 


a : b = r 0 = q l 

(x 5 - 2x 4 + 2x 2 - 1 3x + 1 8) : (x 3 - x 2 + 8x - 20) = x 2 - x - 9 
21x 2 + 39x-162 = r, 


r o : r i = <1 2 

x 20 

(x 3 - X 2 + 8x - 20) : (2 lx 2 + 39x - 162) = 

21 147 

103 Ox 2060 
49 49 ~ Vl 


(2 lx 2 +39x-162): 


103 Ox 2060 


49 


49 


0 = /; 


q 3 

21-49* 81-49 

1030 + 1030 


The last remainder + 0 = r 2 which can be factorized: 


1030 

49 


(jc-2). Hence x - 2 is the 


GCD(a,b). Students can check their operations using DERIVE or TI-NspireCAS as follows: 


# 1 : 

# 2 : 

# 4 : 

#5: 

# 6 : 

#7: 

#S: 


L = 


5 4 2 3 2 

x - 2-x + 2-x - 13-x + IS, b := x - x + S«x - 20_ 


[ql := QUOTIENT (a, b) , rl := REMAINDER^, b)] 

r 2 2 1 

Lql := x - x - 9, rl := 21 -x + 39 -x - 162 J 

[q2 := QUOTIENT (b , rl) , r2 := REMAINDER (b , rl)] 

x 20 1030 -x 2060 

q2 := - , r2 := - 

21 147 49 49 

[q3 := QUOTIENTfrl , r2) . r3 := REMAINDER (rl , r2)] 


1029- x 3969 
q3 := + 


1030 1030 

1030-0 


r3 := 0 

t 


FACTOR (r2) = 


49 


#9: P0LY_GCD(a, b) = x - 2 
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a -x 5 - 2' x 4 +2' x 2 -13' x+18:£:=x 3 -x 2 +8' x-20 

x 3 -x 2 +8' x-20 

q 1 : =po ly Quotient^ b ) 

x 2 -x-9 

r 1 : =poly Remainder^ b) 

21' x 2 +39' x-162 

q 2 : =po ly Quotient.( b,rl) 

x 20 
21 147 

r2 : =po ly Rema inder ( b,rl) 

1030- x 2060 
49 49 

q 3\ =polyQuotient(r 1, r2) 

1029- x 3969 

+ 

1030 1030 

r3\ =po ly Rema inder (r l s r2) 

0 

factor (r2) 

1030- (x-2) 
49 


This cannot be done with the V200 / TI-92. The only one function which seems to apply is 

propFrac: 


It looks fine - at the first glance. The 
quotient of the first division is the 
same as above but the remainder 
does not match. What happened? 

The V200 is too clever! It cancels the 
common factor (x - 2) so it cannot be 
detected any longer. 



A Pretty Knot 


[(3 + SIN(v) + COS(u))-SIN(2-v), 
(3 + SIN(v) + COS(u))COS(2-v), 
SIN(u) + 4-COS(v)] 

0 < u, v < 2n 
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Latest News: 

In the DERIVE-News JiscMail List circulated the following message: 


If we don't have any volunteer then the list will be closed on Friday 26th 
April. When this happens subscribers will not be able to post messages to 
the list but archives will be kept available. 

Thanks for your co-operation. 

JISCMail Helpdesk 


Fortunately a subscriber offered to maintain the list in the future. His first message to the current sub- 
scribers was: 


Dear Derivers, 

My name is Lars Erup; I am a Satellite Communications Systems Engineer based in Mont- 
real, Canada. 

I have agreed to take over the task of acting as Owner for the "Derive-News" list. Let me as- 
sure you that I do not have any plans for changing anything, certainly as long as the list is 
ticking along quietly. Let me also remind you of the current list settings: 

- Anyone can subscribe 

- Only people subscribed to the list can post messages 

- Archives are kept for this list 

- This is a public list, with public archives; anyone can find it (& search archives) on the 
JISCMail website and through search engines 

There are currently 172 of us. 

Best Regards 

Lars 

It s great that this list survives and I find it great that we know the person doing the job. I wrote to 
Lars, offered cooperation and invited him to join the DUG. He answered immediately: 


In your newsletter, you may want to mention the existence of the list (and where to find it), 
and encourage people to use it for discussions. 

Likewise, it might be a good idea if you post a synopsis of the newsletter to the list whenever 
a new edition is ready. 

These are just a couple of ideas. 

Best Regards 
Lars 

The email address of the list is derive-news@jiscmail.ac.uk 

I invite you to use this form of communication between DERIVERs from all over the world 
and subscribe. You can visit the DERIVE-NEWS archive at 

https://www.iiscmail.ac.uk/cgi-bin/webadmin? AO=derive-news 


Call for Submissions 

Deadline: February 15 th , 2014 

Technology and its 
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